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We investigate the generation of primordial magnetic fields from stochastic currents created by 
the cosmological transition from inflation to reheating. We consider N charged scalar fields coupled 
to the electromagnetic field in a curved background and derive self-consistent equations for the 
evolution of the two point functions of the fields, which in the large-N limit give a decoupled set for 
the scalar and the electromagnetic functions. The main contribution to the electric current comes 
from the infrared portion of the spectrum of created particles, and in this limit the damping of the 
magnetic field is not due to normal conductivity but to London currents in the scalar field. For a 
given set of the physical parameters of the problem, we solved this equation numerically and found 
O , that, due to the fact that the London currents are oscillating, the field actually grows exponentially 

• during the time interval in which our large-N limit equations are valid. Although for the chosen 

parameters the induced field is weak, the present uncertainties on their actual values leave open the 
possibility for higher intensities. 

<N 

Cn ■ I- INTRODUCTION 

> , 

One of the most intriguing observational aspects of astrophysics is the detection of large scale magnetic fields in 
\ all the cosmological structures of the Universe . The remarkable aspect of the observations of galactic magnetic 
■ fields, is that their intensity and structure are more or less the same, no matter if the galaxy is one of our neighbors 
or a highly red-shifted one. Clusters of galaxies are also permeated with intense and quite coherent magnetic fields. 

The main puzzle is the dynamical origin of these large scale fields. For galactic fields, two main lines of research are 
pursued: one proposes an origin based on local astrophysical processes, while the other advocates a primordial origin. 
The seducing aspect of a primordial origin of these fields is that it accounts for the observations in all the mentioned 
structures, no matter how high their red-shift is. 

The main difficulty for generating magnetic fields in the Early Universe is the need to break the conformal invariance 
of the Maxwell field. This can be done in many ways, which explains the large variety of proposed mechanisms in the 
literature p|-|2l|l. After this first obstacle is overcome, the success of any given proposal hinges on whether the energy 
from the primordial source may be successfully stored in the magnetic field, or else it is dissipated away (over and 
above the strong reduction brought by cosmic redshift). This point was highlighted in a recent paper by Giovannini 
and Shaposhnikov pl| . Using kinetic theory, these authors showed that the amount of the magnetic field which may 
' be generated during reheating is highly sensitive to the conductivity of the cosmic medium. 

In this paper, we extend this analysis to fields generated in the earliest stages of reheating. At this time the fields 
and their sources span super-horizon scales and so they may hardly be described in terms of more or less localized 
excitations, counted by a distribution function obeying a transport equation. It becomes necessary to recast the 
analysis of p2| in terms of a more fundamental description of the process. We shall atempt to obtain an upper bound 
for the magnetic field which may be generated during reheating, directly from quantum field theory (for a discussion 
of the circumstances in which quantum field theory reduces to kinetic theory see [p3[). 

As a representative mechanism for field generation, we shall adopt the one proposed earlier in ref. J^£4J]. In this 
proposal, conformal symmetry is broken by coupling the Maxwell field to charged scalars fields. The transition from 
inflation to reheating produces a strong amplification of the scalar fields (particle creation) . This results in stochastic 
currents which eventually decay into the magnetic field. The main advantages of this mechanism over others proposed 
in the literature is that it is naturally superhorizoncd and relatively simple, involving no new physics. The bulk of the 
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field is generated early in the reheating era (the so-called preheating period) when on the other hand a quantum field 
theoretic treatment may be atempted. It is realistic in the sense that there are actual candidates for the scalar fields 
in supersymmetric versions of the standard model. We shall ignore the difficulties associated with the fact that the 
electromagnetic field does not exist as such during this period (we should rather follow the evolution of the several 
vector fields in some unified theory, a combination of which becomes the photon after electroweak symmetry breaking) 
as we believe our order of magnitude estimates shall be validated by this more complete treatment. 

Given the present state of knowledge, a quantum field treatment can only be done within some perturbative scheme. 
We shall assume the model contains N identical scalar fields, and thus avails ourselves of the l/N expansion |p5|j26f . 
Working at leading order, we shall obtain the self-consistent evolution equations for mean fields and propagators using 
the closed time path, two particle irreducible effective action (2^J|^] . These tools have been applied to study a broad 
spectrum of non equilibrium problems, from heavy-ion collisions and pair production in strong electric fields p6[ | to 
the study of non equilibrium dynamics of the inflaton field during reheating p8| and to the evolution of quantum 
fields and production of density perturbation in inflationary dynamics [^9| 

We shall study the evolution of a system of N charged, massive scalar fields coupled to one electromagnetic field, 
in an expanding universe which realizes a transition from an exponential expansion (inflationary period) to a less 
rapid expanding epoch (reheating period). For minimal couplin g o f the scalar field to gravity, we can assume that 
during inflation the charged field is in its invariant vacuum state |30|] , there being no particle creation and no electric 
currents. 

When the transition from inflation to reheating takes place, the vacuum state changes and the state of the field 
turns into a multi-particle state. From then on, we have an out of equilibrium system of charged and electromagnetic 
fields, that evolves in an expanding background, and it is this evolution that we shall study. 

The reheating epoch of the universe is very difficult to address, both from the conceptual and from the technical 
point of view. Very little is understood about the process of decay of the inflaton into the other fields to which it is 
coupled, and the consequent establishment of the primordial plasma pl[ | . A complete description of the evolution of 
our system of scalar particles and electromagnetic fields should take into account their possible couplings to the other 
fields present. To atempt to take into account all these effects in a consistent way, makes the problem conceptually 
and technically insolvable. However we can try and incorporate them in a phenomenological way by considering 
a "thermal mass" for the scalar particles, besides their bare one. On the other hand, any other coupling to the 
forming plasma would amount to consider an external electric conductivity, and this case was already considered in 
the literature @§l2l. 



We must observe that to preserve the very weak coupling of the inflaton during inflation 32 1, it may be presumed 
that the inflaton is neutral and uncoupled to other charged species. Therefore the plasma of charged particles would 
not result from the decay of the inflaton, but from some indirect mechanism (e. g., from gravitational particle creation 
p3| ). As a result, the primordial charged plasma contributing to the thermal mass of our scalar fields, need not be 
in equilibrium with the neutral particles, including the inflaton, which contribute most to the energy density of the 
Universe during reheating, and it will take longer to form. As we shall assume that the primordial plasma is not 
established immediately after the end of inflation, we shall consider no external electrical conductivity in the period 
of time during which we follow the evolution of our system of charges and magnetic field: the damping of the fields 
and currents would be due to their own interaction. 

The paper is organized as follows: In the next section we present the general tools that we shall use in the magnetic 
field evaluation: we evaluate in conformal time the scale factors of the Universe during inflation and reheating; we 
write down and solve the Klein Gordon equation for the scalar field modes in each epoch, and calculate the Bogolyubov 
coefficients. As it is not possible find a closed analytical solution to the scalar field equation during the reheating 
epoch, we split the momentum interval in two parts: an infrared and an ultraviolet one and find the solutions for 
each of them. The resulting Bogolyubov coefficients show that the bulk of the particles are created in the infrared 
portion. In section III we write down the Schwinger Dyson equation for the Hadamard two point function of the 
electromagnetic field four potential, D lflv (x, x') = ({A^ (x) , A v (x)}) which will give us the information we seek about 
the induced magnetic field. This is in fact a set of coupled equations for the different components of the two point 
function. They posses two kernels: a local one, and a non local one that can be associated with dissipation. The 
main contribution to these kernels comes from the infrared sector. In this situation, the equations for the transverse 
part of the pure spatial components of the Hadamard function, i.e. the ones whose curl gives the magnetic field, 
decouple. When computing the kernels, the dissipative one turns out to be several orders of magnitude smaller than 
the local one and hence can be ignored. The remaining local kernel will be responsible for screening the field like in 
super-conducting media, i.e. Meissner effect. 

In section IV we translate our equations to a stochastic formulation, which allows a better understanding of the 
field evolution. We write down a Langevin-like equation for the magnetic field itself, and find a numerical solution 
for a given set of physical parameters of the problem. Due to the fact that the local kernel is an oscillating function 
of time we find that, for all astrophysically interesting scales (i.e. the ones that are superhorizoned during reheating), 
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the field intensity grows non adiabatically with the time, a behavior that can be understood in terms of stochastic 
resonance J34|,|35| . Finally in section V we shortly review and discuss our results. 

In Appendix 1 we write down the main ingredients for our study, the Lagrangian for scalar electrodynamics in 
curved space time, we build the 2PI closed time path effective action and from it derive the set of Schwinger-Dyson 
equations for the scalar and electromagnetic two point functions. In taking the large N limit we find that the equation 
for the scalar field propagator decouples from the one for the electromagnetic field, and hence there is no back-reaction 
on the electric current. All the information about the dissipative properties of the system is encoded in the kernels 
of the equations for the electromagnetic two point functions. These equations are a set of coupled equations for the 
different components of the electromagnetic two point function. 

In Appendix 2 we compute the kernels that appear in the Schwinger - Dyson equations for the ultraviolet sector 
of the mode spectrum and for the infrared one, and conclude that the latter dominates over the former. We work in 
natural units, in which c = % = 1, and with signature (—,+,+,+). 

II. SCENARIO FOR MAGNETOGENESIS 

Assume that there exists a charged massive scalar field in the very early universe. Consider also de Sitter inflation, 
and that during that period of the universe the scalar field is in its invariant vacuum state and hence there is no 
particle creation |j0| . When the transition from inflation to reheating takes place, due to the change in the geometry 
of the universe the vacuum state becomes a multiparticle state p7| , which means that a stochastic electric current 
unfolds which induces a magnetic field. 

As stated in the introduction, the reheating period is very difficult to address. It is widely accepted that the general 
picture is the decaying of the inflaton field into other matter fields through non linear oscillations. This is an out 
of equilibrium process, during which the plasma that will determine the evolution of the universe at the subsequent 
epochs is being formed. It is usually assumed that it ends when (conformal) thermal equilibrium is achieved, which 
in turn determines the beginning of the radiation dominated epoch. Very little is understood about the formation of 
the plasma, when the charged matter arises and how the temperature evolves |3~ij ] 

Our created particles will dive through a sea of forming matter fields and it would be very naive to neglect a possible 
interaction between these two systems, no matter how little we know about them. We shall consider that this possible 
interaction produces a shift in the value of the scalar mass, that can be much larger than its bare value. Due to 
the uncertainties about the very process of particle creation (e.g. the change in the geometry is not instantaneous 
and hence neither is particle creation) and about the onset of reheating, we shall take into account this shift when 
calculating particle creation. When studying the propagation of the induced magnetic field, we shall not consider 
a possible conductivity of the forming plasma. On one side this is a reasonable assumption for the early stages 
of reheating, according to the comments made above. On other side this assumption will allow us investigate the 
damping in the magnetic field due to its own sources. 

The scenario we are considering is then as follows: we evaluate particle creation by matching the modes of the 
scalar field and their first time derivatives at the instant of transition between inflation and reheating. We assume 
that this instant is r = and that during inflation the universe expands exponentially. In the early stages of reheating, 
the dominant form of matter is still the oscillating inflaton field, and thus the universe expands like under matter 
dominance |38[|32[| . There being no compelling argument to the contrary, we consider minimal coupling of the matter 
field to the geometry. To model the generation of magnetic fields during this period, we shall follow the interaction 
of the Maxwell field with N identical charged scalar fields, minimally coupled to the Maxwell field. To leading order 
in the l/N expansion, these fields obey the free Klein-Gordon equation (see Appendix 1), and may be decomposed 
into modes in the usual way. We shall assume geometry is described by a spatially flat Friedmann- Robertson- Walker 
model. 



A. Scale Factors of the Universe 

To give a definite form to our model, we must find the scale factors of the Universe for the inflationary and the 
reheating periods in conformal time, which is defined as drj — dt/a(t), t being the physical or cosmological time. 
For the inflationary epoch, the scale factor of the Universe reads aj (i) = aocxp (Hi), so the conformal time is 
r) — rjo = — cxp (—Hi) /Hao and in this variable the scale factor reads ai (i]) = H^ 1 (r}o — f])_^ ■ 

Assuming that during reheating the universe expands as if it were matter dominated p2| , in cosmological time 

the scale factor is or (i) = a\ (t + ti) 2 ^ 3 . Therefore the conformal time reads 77 — 771 — 3 (t + ti) 1 ^ 3 ja\ and the 
corresponding scale factor an (rj) = (77 — 771) /9. 
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In order to avoid overproduction of created particles, when matching the scale factors we must demand that at 
the time of transition, the scale factor changes smoothly. The minimum conditions to be fulfilled are that the 
functions and their first derivatives be continuous at the transition time, then obtaining ai (77) = (1 — Hrj) 1 and 
an (77) = (1 + if 77/2) . Defining the dimensionless time r = Hr) the corresponding functions read ai (r) = (1 — r) 
and a R (r) = (1 + t/2) 2 . 



B. Field Equations and Bogolyubov coefficients. 

In this section we solve the Klein Gordon equation for the scalar field, in the two eras of the Universe involved 
in our study and calculate the Bogolyubov coefficients. From now on we will work with dimensionless variables 
k = Hk, m — > m/H, together with the dimensionless time variable r defined above. 

Assuming minimal coupling to curvature the Klein Gordon equation for the modes of the scalar field reads 



<r , 2 m 2 (t) 2 a (t) 

+ k + —ET 1 ^ ( r ) - 



fk (r) = (1) 



dr 2 H 2 w o(t)_ 

In spite of our ignorance about the reheating epoch and the actual interactions of our scalar field during it, it would 
be very naive to neglect them, so as stated in the Introduction, we shall take them into account by considering a 
thermal mass given by the phenomenological expression 

m 2 (r) = gT 2 (r) (2) 

with g being a coupling constant which we consider of order one. We assume a generic form for the evolution of the 
thermal mass during reheating, i.e. one of the form 

a b (r)" 

where 6 is a parameter that satisfies < b < 1, and Tm is the maximum value of the shift, attained at r ~ 0. This form 
of the mass shift is inspired in the process of preheating, where a peak temperature is attained which subsequently 
decreases as reheating unfolds. Setting b = 0, g = 1 implies Tm = 7no, where mo is the bare mass of the scalar field. 
For any other values of b > it must be Tm > ma- 

The Bogolyubov coefficients are given by the decomposition 

fk(T)=a k f?(T)+f} k f?(T) (4) 

with (r) the modes of the scalar field during reheating, and can be obtained by demanding continuity of the mode 
functions and their derivatives at the time of transition. 



r(r) = ^ ( 3 ) 



1. Inflation 

During this period, as the temperature of the Universe is practically zero, the modes satisfy the bare mass equation, 



ml/H 2 



dr 2 



(i- T y 



fk (r) = 



1/2 

Writing f k (t) = (1 — r) h k (r) we obtain a Bessel equation for h k (t) 



1 



dr 2 



(1 - t) dr 



+ k 2 - 



9/4 -m 2 /H 2 



h k (t) = 



(5) 



(6) 



whose positive frequency solutions are the Hankel functions of the first kind, [k (1 — t)] with v = 

(3/2) \/l — 4mg/9iJ 2 ~ 3/2; the last relation stems from the fact that during inflation rn^/H 2 <C 1. We then 
have that the normalized, positive frequency modes during inflation read 



(l-r) 1/2 H^ 2 [k(l 



,ik(l-r) 



2k 



fc(l-r) 



(7) 



The last equality follows from the exact polynomial expression for H^ 2 [k (1 — r)] |3t 
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2. Reheating 



For this epoch the Klein Gordon equation reads 



d T 2 



T 2 
H 2 



k 2 + g±K [ 1 



4(1-6) 



2(1 + 7-/2)' 



A (r) = 



(8) 



In this case it has no closed analytic solution, so we shall solve it in two limits, namely the large wavenumber 
limit, for which k 3> A = g 1 ' 2 Tu/H and the small wavenumber limit, for which k <C A. Care must be taken when 
doing this, because in eq. (|) the factor g^ 2 T M /H is multiplied by a growing functions of time. This means that in 
principle this splitting is time dependent. However this will prove to be unimportant at the end of the calculations. 

a. Large wavenumber limit The Klein Gordon equation in this limit reads 



JUL. 

dT 2 



+ k 2 - 



2(l + r/2)' 



(9) 



proposing again f& k (r) = (1 + r/2) 1 ^ 2 h R ^ k (r) we obtain a Bessel equation for hh (r): 



h R(l)k (r) 



1 



(l+r/2) 



fiR(l)k ( T ) 



Ak 2 



9/4 



h R {i) k (t) = 



(10) 



( 2) 

whose positive frequency solutions are the Hankel functions of the second kind [2k (1 + r/2)J. The normalized 
mode functions read 



/&M=yf( 1 + I)'% : 



(2) 
3/2 



2fc 1 



3 -i2fc(l+r/2) 



1 - 



2fc(l + t/2) 



(11) 



where in the last line we have used the polynomial expression of the Hankel function. 

Replacing the corresponding modes and their derivatives in the general expression for the Bogolyubov coefficients 
we get 



Oi(l)k 



Slk 
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3/ 



— 

+ 2k 8fc 2 8fc 3 



(12) 



0(i)k = e 



— ik 



3 _3i_ 

8k 2 ~ 



(13) 



b. Small wavenumber limit For small wavenumbers, i.e. those for which k <C g 1 t 2 TM / H , the equation reads 



M 



dr 2 ' 9 lP l 1+ 2 



T \ 4-46 



1/2 



/( S )fc (t) = 



Assuming as time variable u = (1 + r/2) and replacing (z) = it 1 / 2 /! (w) we get 



(14) 



1 ■ 



h( s )k (u) + -h {s)k (u) + 



T 2 
H 2 



h (s)k (u) = 



with c = 4 — 46. Defining x — u 1 the equation now reads 



^ 2 , / \ 1 d i , \ 



4.9 



#2 7 2 



47 2 x 2 



(15) 



(16) 



which can be cast in the form of a Bessel equation by demanding c/7 + 2/7— 2 = 0. We get 7 = (3 — 2b) and the 
equation reads 
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-^ h is)k{x) + - — h {s)k { X ) + 



4ff 



T 2 1 

1 M 1 



H 2 7 2 47 2 x 2 



h( s )k (x) = 



(17) 



whose positive frequency solutions are again Hankcl functions of the second kind, Hf/ 2i [(2g^ 2 T M /Hj) u 1 } . The 
normalized solutions to the field equation are then 



with 



z(t) = z (l 



T\7 



Z = 



23 1 / 2 T, 



M 



2J ' j H 

Replacing in the general expressions we get by a straightforward calculation 



7r e 



ik 



27\/2fc 



1 + i 



n l/2 T M „(1) r i ,-1.17(1) 1 



3/2 7 I 



(18) 



(19) 



(20) 



Ac = -i 



7r e 



27 V2fc 



1 + 



l/2^M ^(2) r 1 

9 ' -Jf H \/2 1+ i N 



and in the limit zo « 1, which is always valid, they reduce to 



3/27 



ik 



k 3/2 



(21) 



(22) 



Observe that the apparent breakdown of these formulas at 7 = (due to the inadequacy of the Bessel function 
representation of the solution) lies outside the physical range. 



III. DYNAMICS OF THE ELECTROMAGNETIC TWO POINT FUNCTIONS 

Cosmological particle creation is a stochastic process. If the particles are charged, electromagnetic fields are induced 
by the created currents. But their mean values being zero, they manifest through their variances or two point functions. 
The functions we need to obtain information about the evolution and the state of the field are respectively the retarded 
D r J* (x,x') and the Hadamard D\ V1 (x,x') two point functions, defined as 



D r * (x, x 1 ) = 1 ([A v (x) , A 1 (x')]) 6 (r - r') 



(23) 



D luj (x, x') = ({A v (x) , A 1 (a/)}) 
The evolution equations for this propagators, known as Schwingcr-Dyson equations (see Appendix 1), are 

»T + (l - J) d»8 v - e 2 T^ (x, x) D% (x, x') 



+^e 2 J dx" (x, x" ) D r u e j (x" , x') = - 5^5 {x - x') 



rT^x ^ - e2r ii ( x < x ) 

+ie 2 J dx" £{£ t (x, x" ) D lvy (x" ,x') = - 
where ( is a gauge fixing constant, U x = — d 2 + V 2 , and 



-D11/7 (x,x') 
2 J rfx"Sr (x,a?)D% v (x",x') 



e 
~2 



(24) 



(25) 



(26) 



G 



r%(x,J)=Tr [G\ d (x,x)+G 2 od (x,x)] 



(27) 



V£ M = V m {x)^(x") [G\ c , (x,x»)dJpG dd , (x,a»j\ (28) 

c and d are closed time path indices whose values are 1 for the forward directed time path and 2 for the backward di- 
rected time path and the supraindex i = 1, 2 denotes the real and imaginary parts of the complex scalar field. According 
to the values of the closed time path indices we have that G21 (x, x 1 ) = {<$> (x) 4> (%')} is the positive frequency two point 
function for the scalar field and G12 (x, x') = (<f> (x 1 ) <fi (x)) the negative frequency one. With these basic propagators we 
can build the antisymmetric two point function, also known as Jordan propagator, G (x, x') = G21 (x, x') — G12 (x, x') 
and the symmetric, or Hadamard one, Gi(x,x') = G2i(x,x') + G12 (x, x'). Gn(x 1 x') = (T (<f> (x) (p (x'))) is the 

Feynman propagator and G22 = \T {(f> (x) <t> { x '))) the Dyson one. The retarded and advanced two point 

functions are defined as G re t(x,x') = iG[x,x') = G a d v (x',x), or else G re t{x,x') — i [G11 (x, x') — G12 (x, x')], 
G a dv {x,x') = i [G22 {x, x') — G\2 (x, x')]. The same definitions apply for the two point functions for the electromag- 
netic field, i.e. D 2 ^ (x, x') = (A^ (x) A v (x')) is the positive frequency two point function, D 12 (x, x') = (A^ (x') A v (x)) 
the negative frequency two point function, and so on. Based on these definitions, we have built the non local 
1) and (g), as Ej£ (x,x») = S^,n (x,x») - S^ 12 (x,x») = S^, 21 (x,x») - S^ 22 (x,x") and 



kernels in eqs. 

^1 (x, x ) = 5D 1212 {x> x ) + E 2121 [x, x ) = E^ 1;11 (x, x ) + S 2222 [x, x ). 

Each of equations ( p5| ) and (]2^) is in fact a set of coupled differential equations for the different components of the 
retarded and Hadamard two point functions of the electromagnetic four potential. Equation (^6|) is the one that we 
shall use to evaluate the magnetic field. 



A. Equation for the spatial components of the electromagnetic two point functions. 



The equations for the electromagnetic two point functions form a set, where the equations for the pure spatial and 
pure temporal propagators are coupled by the mixed functions, i.e. by two point functions with one temporal and 
one spatial component. This coupling is realized through the non local kernels, specifically by its — i components, 
which in general do not vanish. 

When we take into account particle creation, we see that local and non local kernels split into a vacuum and a 



particle contribution, i.e. (x, x) — (x, x) +P 



(x,x) and £^(1) (x, x" 



We are interested in the contribution of the created particles and hence shall solve 



- l ret(l) 



(x,x»). 



1 



( 1 - ^ ) o^r - e 2 r^ (/,; ( 3 .,s 



D Xvn (x,x') 



(29) 



-ie 2 / dx" E^ (P) (x, x" ) D lvi (x" , x') = - & - / dx" ^ [P) (x, x" ) (x" , x') 



For our purposes of magnetic field evaluation, we shall need only the spatial components of the propagator 
Duj (x,x'), i.e. the solution to the i — j equations, and in particular its transverse parts. Since the bulk of particle 
creation occurs for long wavelengths, the main contribution to the kernels comes from that sector of the spectrum 
(see Appendix 2). The — i components of the non local kernels take the form of a gradient and in that limit the 
mode functions depend on the momenta only through their moduli. When transforming Fourier the spatial part of 
the two point functions, there remains an integral over the momenta that vanishes for the transverse part of the — i 
component. This means that the equations for the transverse part of Duj (x,x') decouple from the set, a fact that 
facilitates enormously their resolution. The Fourier transformed equations that we shall work with then read (see eq. 
(p6|) and Appendix 1) 



n a (d 2 + k 2 )+e 2 



d 3 p 
(2^) 



3/2 



(30) 



-^75 / dr" (p, k-p,r, r" ) D UJ (fc, r" , r') 

(2tt) j/ J 

y J dr"Ef P) (p,k- P ,r,^)Dtf v {k,r',r") 
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where k is the comoving wavenumber of the spatially Fourier transformed two point function of the electromagnetic 
field, and p the corresponding one of the scalar field. For example when replacing the mode decomposition ( |4j) the 
local kernel reads 

Tn (P) (P,r,r) = rf l {2 \(3 p \ 2 \f p (r)| 2 + a p (3* p f 2 (r) + (3 p a* p p p 2 (r)} (31) 

with f p (r) given by eq. ( |Ts| ) or ( |ll"|) and the Bogolyubov coefficients by (|2^), and similarly the non local kernels. We 
assume that the vacuum part may be absorbed into a renormalization of the classical action, the remainder being 
negligible. 



B. Computing the kernels 



In this section we shall compute the kernels found in the equations (|3C|) for the Maxwell field propagators. It may 
be checked that the contribution from the large wavenumber sector is negligible (see Appendix 2) in comparison to the 
infrared one, a fact that can be understood by looking at the expressions for the Bogolyubov coefficients, eqs. jl2]), 
(H) and (|§). The mode functions are given by eq. ( J18| ) an d we see that they do not depend on the wavenumbers 
p. The Bogolyubov coefficient on the other side retain their full momentum dependence. 



When we replace for the modes and coefficients, eqs. (18) and (|22[), we find a logarithmic divergence in the 

momentum integral of the local kernel eq. (|3l|). We introduce an infrared cut-off, which we can choose as the 
comoving, dimensionless mode T corresponding to the original inflationary patch that gave rise to our universe. For 
the upper limit we take A ~ Tm/H. We then have 

r n t) = r (p. t > - ^ . ( 32 ) 

Jr (2n) i/2 

f = 2 3 ^r 2 (*L±3) In (£) z^F 2 (z) (33) 

where 

F(z) = z 1 ^J 3/2j (z) (34) 



J3/27 (z) being a Bessel function of the first kind an z and zq being defined above in eq. (|19[). 

The non local noise kernel has no infrared divergences and is A 2 times smaller than the local one, with A ~ 
Tm/H <C 1 (see Appendix 2). This difference in the orders of magnitude will never decrease, (even if we evaluate 
the remaining time integral in the non local part of the equation) and therefore for our purposes we can ignore the 
contribution of this kernel. 

This is a very important approximation, not only from the technical point of view but also because it shows that 
dissipation due to the ordinary conductivity of the own charges is negligible. 

We now turn to the Hadamard noise kernel T,[ j(N) (p, k — <p,T,T n ) that is the source in the equation (^3) for 
d[ N) (k, T ,T'). The same considerations made for the retarded kernel about the gradient structure of the — i 
component of the two point function hold for this one. So we shall only need to evaluate the pure spatial components, 
given by 

£f P) (p,fc-p,r,T") ~ PV |3 ft fr,^), (35) 
p 3 \ k — p 

Si (r, t") = 2 6 /^ 1 r 4 (p^j *o Sh F 2 [z (r)] F 2 [z (r")] (36) 

It is remarkable that the same function F appears in this kernel and in (|33|). 

We see that after replacing the modes and Bogolyubov coefficients the resulting kernels given by eq. ( p3] ) and ( |3"6| ) 
are real functions. This happens because of exp. (^2|) for the Bogolyubov coefficients and the fact that the modes 
( |l8| ) do not depend on the wavenumbers. In other words, after replacing the Bogolyubov coefficients, the modes ( |l8| ) 
combine to give real functions that oscillate coherently, an indication that they are indeed superhorizoned and hence 
frozen. 
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C. Evolution equations for the magnetic field 



After all the considerations we have made, the equation for the transverse part of the pure spatial Hadamard two 
point function for the electromagnetic four potential reads 



rf l (d 2 + k 2 ) + e 2 lf i (F) (r, t) D Uj (k, r, r') = Sy (k, r, r') 



with 



E ij (k,r,r') = ^ I dr" I J^Zf( P ,k-p,T,T»)Dlf(k,T',T») 



2 Jo Jt (2tt) 3/2 
Let us introduce the Hadamard two point function for the magnetic field: 



({B t (T,r),B 3 (r',f')}) = H^ kl e^ 1 ' 



/d 3 k 
^-prhk^Diw (k, t, t') exp [ik. (f - f')l 
(2nr /2 



(37) 



(38) 



(39) 



where the prime refers to the prime coordinate and iJ 4 gives the dimensions. 

We are interested in the field Bi (r, A) coherent over a scale A, so we must filter the high frequency contribution 
with a window function of "size" A, i.e. we must evaluate 



1 



d 3 r / d i r'({B i {r,f),B j {r',f)}) 



where V\ ~ A 3 is a comoving volume in which we seek homogeneity. Eq. (|40|) amounts to calculate 



^- [ cfVexp [ik.f] = — / d 3 r' exp [-ik.r ] = W\ (fc) 
Va Jvx V\ J Vx 



This window function can be approximated by 



W x (k) ~ 1 i/ k<K = 1/A 
= otherwise 



(40) 



(41) 



(42) 



which can be implemented as a cut-off in the fc-momentum integral. 

For our purposes it is enough to compute the self-correlation ({B l (t, A) , Bi (r', A)}) . From the equation for the 
Hadamard propagator we get 



^r(^) 



({& (t,X),B 1 (t\X)}) = 



(43) 



2e 2 H 4 



K 



d 3 k 



o (2ttH z Jt (2tt) 



/" j // | fc X | a I in n\ r^ret / I ii\ 

dr — —Si (r , r )D (t , r ) 



,3/2 



p 3 |fc — p\ 



where we used the fact that, because K -C 1, the spatial gradients are negligible, and D ret is therefore independent 
of the wavenumber. The momentum integral can be easily evaluated and we obtain 



dr 



2 +e 2 T(r,T) 



{{B l (r,A),A (r',A)}) = 2e 2 tf 4 A 4 / (r,r")^ ret (r',r") 

Jo 



(44) 



IV. EQUIVALENT STOCHASTIC PROBLEM 

Rather than solving directly equation (|44|), we move to a physically more transparent language and translate our 
problem to a stochastic formulation. We introduce a stochastic field B (r) obeying the Langevin equation 



dr 2 



+ e 2 r (t, r) 



B(r)=f(r) 



(45) 
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where £ (r) is Gaussian noise with zero mean and self correlation 

^(rJ^/J^^eWS! (t,t") 

It is then easy to see that 

({^(r.A),^ (r',A)}> = ( J B(r)A(r')) € 
Moreover, because of the structure of the self-correlation, we may realize the noise as 



f (t) = eH 2 K 2 2 3h T 2 



(46) 
(47) 

(48) 



where X is a single gaussian variable with (X) = and (X 2 } = 1. 
Since the equation is linear, we may write 



B (r) = etf 2 X 2 2 3 ^L 2 (^^) ^ 4/7 *5 s (r) 



The equation to be solved is then 



where 



B s (r) = F 2 [z(r) 



C 2 = i 2 z- 4h 



(49) 



(50) 



(51) 



where e 2 2 3 ^- 1 r 2 [(27 + 3) /2'y] In (A/T) = e 2 with the boundary conditions B s (0) = S s (0) = 0. Eq. @ resembles 
the London equation in a superconducting medium pQ] . This means that the induced field will be weaker than the 
obtained in the case of free propagation. We may write B s (r) = C -2 (1 + Bh (r)) and the equation for i?^ (r) reads 



dr 2 



+ C 2 F 2 [z (r)] 



B h (r) = 



(52) 



with initial conditions Bh(0) = — 1, .B^ (0) = 0. It is convenient to adopt the variable z as independent variable. 
Equation (K2) now reads 



| S 2 -6/7 /„,2 



" d 2 


+ - 


cfe 2 


= 2T 


2(7 



7Z dz 



+ c 2 g 2 (z) 



B h (z) = 



(53) 



given by 



B(r)~H 2 K 2 (l + B h (T)) 



(54) 



Because of the oscillatory behavior of g 2 (z) in equation ( |53| ) we expect that the attenuation of the induced magnetic 
field due to the London currents will not be so strong. Besides, due to the form of g 2 (z) we might expect some sort 
of resonant growth. This can be understood as follows. 

Defining 



and replacing in eq. (j53J) we get 



b h (z) + 



B h {z) = z-^- x ^b h {z) 



7 2 -l 4S 2 Jf /27 (*) 



bh (z) = 



4 7 2 Z 2 7 2 z 6/7 z (27-3)/7 

As in general it is satisfied that 4e 2 z 6 ^ 7 ^> (7 2 — l) /4 we can approximate equation ( f56] ) by 



(55) 



(56) 
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b h (z) + c 2 ^l^ b h (z) = (57) 



Consider the asymptotic form of the Bessel functions, for z ^> 1, i.e. J v (z) <~ \J2/ (ttz) cos [z — ^7r/2 — 7t/4]. The 
second term in eq. (|57j) can be then written as q (7, Tm/H, z) cos 2 [z — (3/27) n/2 ~ it /4] with 

422 

q(y,T M /H,z)~— ^j-z^h (58) 



For g (z) ^> 1, or equivalently until 



4c 



,-2 



2 6 /7 



7/(37~3) 

) (59) 



the growth of the field intensity could be described in terms of stochastic resonance (originally developed in Ref. |34| ] 
and extended in Ref. J35|), because of violation of the adiabatic condition. 

In order to obtain an accurate estimation of the magnetic field, we shall integrate eq. ([53to numerically and evaluate 



the induced field from eq. (54). This procedure will also allow us to confirm the previous heuristic analysis of 
stochastic resonance. In view of the uncertainties in the actual values of 7 and Tm we shall not atempt to estimate 
the induced field for a wide interval of their possible values, but shall illustrate the effect by considering one specific 
choice, namely 7 = 5/2 and Tm/H = 10~ 2 . We also fix e 2 = 1 for all possible values of 7 (in fact e 2 is a slowly varying 
function of 7 and this choice does not affect the outcomes of the numerical integration) . This value of 7 corresponds 
to b = 1/4 and describes a strong deviation of the thermal mass from conformal thermal evolution. We recall that 
the parameter b does not describe the evolution of the temperature of the reheating plasma, i.e. the one formed by 
the inflaton decays, but the one of the thermal shift of the scalar charges. For this values of the physical parameters, 
according to eq. ( |59"| ) we might expect resonant growth until z « 500. 

Considering the maximum allowed value for the Hubble constant during inflation, namely H = 10 13 GeV, Tm/H = 
IQ- 2 corresponds to T M = 10 11 GeV. 



A. Limits of validity 



Eq. ( p3| ) correctly describes the evolution of Bh during a time interval in which the equations obtained in the 
leading order of the 1/N approximation, eqs. (^5|) and ( |26| ) are valid. To estimate this interval we can compare 
the growth of the electromagnetic energy density with the one of the scalar field: assuming efficient conversion of 
scalar energy into electromagnetic one, the moment at which the latter overtakes the former can be considered as the 
limiting time to integrate the equations. In fact, in the leading order of the 1/N approximation the equation for the 
scalar field propagator shows no back-reaction from the magnetic field (see Appendix 1). This means that to this 
order the induced field is very weak and hence its energy density must be smaller than the one of the scalar field. 
When this condition is broken, the evolution equations for the two point functions cease to be valid. The energy 
density is given by the 00 component of the stress energy tensor of the system of fields, T 00 . To evaluate this tensor 
component we consider the electromagnetic field as classical and the scalar one as quantum. We can split T 00 in a 
pure scalar contribution and an electromagnetic part as 



^00 

v 



H 



2 — (pip - 
a 



1 



rf diipdjip ■ 



m 2 (t) 

~1p~ 



a 2 (r) + - 1 * - 



(60) 



rpOO 



H 4 



-H* 



F ij Fn + 2F 0l F< 



(V 



(61) 



rf j eAi (djpip 2 ~ fidjp 2 ) + -rf J e 2 AiA 3 (p\ + tp 2 2 ) 



-"OCA. 



where it is understood that all quantities between brackets are dimensionless. We shall compare \T®°) with {T\~ j: 
the instant at which the the latter equals the former will be considered as the time at which the equations obtained 
in the large limit cease to be valid. 
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1. Stress energy tensor for the scalar field 



Considering the electromagnetic field as classical and the scalar as quantum, and defining T 00 = H 4 T 00 we can 
write 



rpOO 



d 3 p 



jj5 \ <Pp ( T ) ¥>p ( t ) - - [<P* P 0) <i>v (t) + <P* P (r) <P P (r)] 



(62) 



m 2 (t) d 2 



+ r] l3 p t pjif* p (r) ^ p (r) + ( a g2 + ^ ) M (t) 



Writing the modes as cp p (r) ~ a p f p (r) + /3 p /p (r) with the modes f p (r) given by eq. fllS| ) and the Bogolyubov 



coefficients by eq. A22J), the contribution from the created particles to (T^° ) reads 



m00\ 



A j3 



/9 Jt (2tt) 



3/2 



i/9pi 2 /w+rw 



2- 

a 



/W + r(r) [/(r) + .r(r)] 



(63) 



m 2 (t) d z 



>7>^ [/ (r) + /* to]' + ( + ^ ) [/ (r) + /* (r)] 



Neglecting the gradient term because it will give a negligible contribution in comparison to the other terms and 
performing the momentum integral, we get 



rpOO\ 



4?r 7 -3/ 7 v , 

z n x < a 



(64) 



/(r)+r(r) --[/( T )+r(r) 



where we have used a 2 m 2 (r) /i/ 2 = o 1 ~ 1 Tm / H . 

Replacing the modes in terms of the Bessel functions through eq. (|lj), the scale factor and its derivative in terms 
of z and zq and using Tm/H = jzq/2 we are left with 



^ 00 ) - tSt^^t? { j 3/2t [* W] + 4 27+ i [* (-)]} 



(2^)- 



(65) 



For z 3> 1 we can replace the asymptotic expressions of the Bessel functions p9[| , i.e. [z (r) 

\/2/ (7rz) cos [z — ^7r/2 — 7r/4] and so obtain the final expression for IT® ): 



T 



mi 



2 1-1/7 



47 2 z 



,3/2 ,2/7 



(27T 



(66) 



which is plotted in Fig. 3. 



2. Stress-energy tensor of the electromagnetic field 

As stated in eq. (|6l|), the stress energy tensor for the electromagnetic field consists of the pure electromagnetic and 
the interaction energy density. Proceding with eq. ( |6l|) similarly as for T 00 we are left with 

n°) = l(B 2 + E 2 )+e^ Z -^Jl /27 (z) (67) 

Using that B = B/H 2 and 



~ _ d B _ 7 zl h dB 

E ~ ~faK ~ 2K^h^ ~dz~ (68) 
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we get 



t5° ) = ~ 



D 1 



7 2 4 H (dB 
4K 2 z 2 h~ 2 \ dz 



,B 2 z l h 



— —J 2 (A 



(69) 



We solved numerically eq. ([53|) and with the outcomes reconstruct B from eq. (|54| ) . Eq. ([54|) gives us the field coherent 
over a comoving scale K' 1 in which we might be interested. We see from that equation that the field intensity grows 

as K 2 . In order to make a honest computation of (t^ 1 ^ we must take into account the highest possible magnetic 

field intensity produced by our mechanism, and this is produced at a scale if max ~ Tm/H, i.e. the cut-off that we 
used to split the momentum spectrum. In Fig. 3 we have also plotted the outcomes of these calculations. 



V. MAGNETIC FIELD EVALUATION 



To estimate the intensity of the induced magnetic field, coherent over a physical scale k 1 at a given time r, we 
must replace 



a (t) k = HK (70) 

So, from eqs. (JsiJ) we get 

B phys (r) = ^ ~ K 2 (1 + B h (r)) (71) 

With the parameters given above, i.e. 7 = 5/2 and Tm/H = 10~ 2 and fixing e 2 = 1 for all values of 7, we have 
numerically solved eq. (|53|) and show the results in Fig. 1. We see that Bh [(z)] oscillates with an amplitude that 
increases exponentially. We have estimated numerically the exponent of the exponential envelope in the z interval 
considered, finding /i = 2/15. The oscillatory shape of the curve clearly suggests that the growth of the field is 
non-adiabatic, and hence we can expect stochastic amplification for all values of the physical parameters and during 
a time interval such that q (7, Tm/H, z) >• 1. 



4000 - 



2000 - 



-2000 - 




100 



FIG. 1. Plot of Bh vs. z: the magnetic field intensity grows exponentially, as is shown by the dashed line. The numerically 
estimated exponent is fi = 2/15 and the amplitude Boh = 300 



The fact that the field grows non adiabatically can be seen in Fig. 2, where we have plotted B^ and g(z) vs z. We 
can clearly seen that the maximum growth of the field occurs when g(z) = 0. 
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FIG. 2. Plot of Bh and g(z) vs. z: The dashed line corresponds to Bh, as in Fig. 1, while the filled one to g(z). We see that 
Bh grows appreciably when g(z) = 0, a fact that confirms the non adiabaticity of the induction process. 



To calculate we considered the exponential envelope of the numerical evaluation and its derivative to estimate 

the magnetic and electric fields respectively. It can be seen from Fig. 3 that at z ~ 95 the magnetic energy density 
equals the one of the scalar field and hence our equations cease to be valid. At this time, (1 + Bh (r)) ~ 10 3 . The 
growth of the field can therefore be understood in terms of stochastic resonance since, as stated in the previous 
subsection, we can expect this behavior for zq < z < 500, and the allowed time interval of integration of our equations 
is well inside this one. 
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FIG. 3. Plot of the energy density of the scalar field (solid line) and of the electromagnetic field (dashed, oscillatory curve) 
as a function of z. The magnetic energy overtakes the scalar one at z ~ 95 



For example, if we consider the comoving galactic scale today, k ~ 10 GeV and use the equivalence lGeV 

jgal 



10 



211 



Gauss, we obtain a magnetic field intensity (T to d) — 10 Gauss, a value too weak to seed the galactic dynamo 



VI. DISCUSSION 



In this paper we have evaluated in a self-consistent manner the induction of magnetic fields from electric currents 
created at the inflation-reheating transition. As the mean value of the current vanishes, this current arises due to the 
quantum and stochastic fluctuations around the mean. 

Using techniques specially suited for studying quantum fields out of equilibrium such as the 2PI closed time path 
effective action together with the large N approximation, we have obtained consistent evolution equations for the two 
point functions for the scalar and electromagnetic fields, known as Schwinger - Dyson equations. 

To evaluate the field we took into account the possible couplings of the created particles to the other fields present in 
the reheating epoch, by using a phcnomcnological expression of a temperature dependent mass. Any other coupling to 
the forming matter fields could be account as an external conductivity, which would give ordinary dissipation already 
discussed in the literature jl5],|2|j2|] . 
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The first important result of our work is that at leading order in the 1/N expansion there is no back-reaction of the 
induced magnetic field on its sources, and hence the charged field propagates as a free field. 

When computing particle creation, we find that the largest number is created in the infrared portion of the spectrum, 
giving this sector the main contribution to the different terms that enter in the Schwinger Dyson equations. 

The evolution equations for the magnetic field two point function show two kernels, a local and a non local one. 
When we evaluate these kernels with the infrared modes and Bogolyubov coefficients, we find that the local kernel 
dominates over the non local, dissipative one, by several orders of magnitude during all reheating. Physically this 
means that dissipation in our system of fields will not be due to an ordinary electric conductivity. This fact also 
allows the equations for the pure spatial components of the electromagnetic two point functions to decouple, a fact 
that facilitates enormously their resolution. 

In order to estimate the induced field we translate our equations to a more familiar and physically clearer Langevin 
equation, which due to the presence of the local kernel, resembles the London equation for a superconducting medium. 
This means that there exists attenuation produced by the same created particles, which is not due to ordinary 
conductivity but to screening due to the fact that the created particles behave as a superconducting medium. This 
is the second important result of our paper and could be obtained because we work with the full quantum field 
theoretical features of our system. 

To visualize the process of field induction, we have solved numerically eq. (|53j) for 7 = 5/2 and Tm/H = 10~ 2 
and found that, due to the fact that the London current oscillates, the damping is not perfect. The magnetic field 
grows exponentially, and the growth rate in the analyzed case is /i ~ 2/15. This exponential growth is non adiabatic 
and generic for all superhorizoned modes of the magnetic field, and could be cast in the framework of stochastic 
resonance. In view of our ignorance about the actual couplings of our scalar field and specially about the very process 
of reheating, we did not pursue a thorough analysis of the possible parameter space in order to give bounds to the 
final field intensity. 

In general we can estimate the upper limit of integration of eq. (^3|) as the one at which the energy density in the 
electromagnetic field overtakes the one in the scalar field. This estimation is based on the fact that to leading order 
in the 1/N expansion, there is no back-reaction of the electromagnetic fields on their sources. When this condition is 
violated, i.e. when the electromagnetic field becomes so intense that its energy density equals the one of its sources, 
the approximation breaks down and the equations cease to be valid. 

For the set of parameters used in this paper to illustrate the field induction, this time interval is rather short and 
the resulting field intensity too weak to have an astrophysical effect. 

Due to the different powers of the momenta, the kinetic term in eq. (^) for the 00 component of the electromagnetic 
stress-energy tensor is negligible with respect to the interaction term, i.e. e 2 A 2 (ip 2 ). The growth of this term is given 
by the growth in the magnetic field, which in turn is determined by the index of the exponential. Therefore the 
smaller this index the later the magnetic energy will catch up with the one of the scalar field, which means that the 
final field intensity will be higher. 

D. Boyanovsky et al. have shown in Ref. |^l| that the electrical conductivity of a primordial plasma might grow 
exponentially. Therefore it could happen that outside the range of validity of the 1 /N approximation, the conductivity 
might damp the amplification of the field. The final answer, however, can only be given from the analysis of the higher 
order corrections to the magnetic field evolution. 

Given our present ignorance about the actual values of the physical parameters that determine that index, a very 
small index and hence a stronger field cannot be ruled out. Also the obtained intensities could be preamplified by 
some primordial process, as for example the one considered in Ref. j4^j4j|, and then lead to acceptable values to seed 
the subsequent amplifying mechanisms. 
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VII. APPENDIX 1: SELF CONSISTENT SCALAR ELECTRODYNAMICS: CLOSED TIME PATH TWO 
PARTICLE IRREDUCIBLE EFFECTIVE ACTION AND 1/N APPROXIMATION 

In this appendix we derive the self-consistent, causal equations for the scalar electrodynamics, in the large N 
approximation, using the closed time path formulation of quantum field theory in a Friedmann Robertson Walker 
universe |p7| , p8| . This formulation helps to correctly identify the Feynmann graphs that contribute to each order in 
the 1/N expansion of the equations. 
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The starting point for its derivation is the action for a set of N charged scalar field coupled to the electromagnetic 
field, which in curved spacetime reads 

S F 4,A^j = - J d^x^f—g {g»» - ieAJ fa (d v + ieA v ) <j>\ + (to 2 + X R 2 ) <t>i<t>\ (72) 

where x = (t, f), e is the electric charge, g^ v the metric tensor, y/—g the determinant of the spacetime metric, \ t ne 
coupling constant of the scalar field to gravity and £ a gauge fixing constant. Sum over repeated greek indices is 
assumed. The complex field can be expressed in terms of real ones as <f) — [<f>i + i<h] I V% and ft — [4>\ — ifo] /V%- 
To implement the large N limit it is convenient to re-scale the fields and coupling constant as 

ti^VNti, — > y/NAp, e^-^= (73) 

obtaining 

S^ = -Nj d 4 xV=9 { Ig^d^i&ti + \ (to 2 + X R 2 ) (<t>? + (74) 

+ \g» a g vP F aP F^ - ±g^g^d^A a d„A f ^ 
which means that each vertex appears N times in the interaction. 



A. Closed time path generating functional for N scalar fields in curved spacetime 



The closed time path generating function for the n-point functions is obtained by coupling the scalar and electro- 
magnetic fields to sources J a (x) in each of the time paths 

r N 

Z [J, K]=Z G \{ DA%_D(j>i±D<f)%± exp {iSg [0?±, <^ ± , 

J 71=1 

+z J d^x^—gC ab [jf (x) C (x) + J£ (x) A, b (x)] 
+ \J d A x^g~ J d 4 x'y/^g 1 C ab C de 

x ~KZ m x') 4>T (x) 4P e m (x) + K% (x, x>) A, b (x) A ve (x') 



with 



+ d 4 x^J—g I d A x'^C ah C d& 

x KlJ (x, x 1 ) [C (x) A, e (x') + C (x') A, b (x')] } 
exp [iW (J, K)] 



s" [m ± ,m ± ,a» ± ] = S" [fi + ,<tq + ,A$]-S» [0?_,^_,< 



(75) 



(76) 



where the plus (minus) sign corresponds to forward (backward) time direction, n,m — 1 . . . N. and i, j = 1 or 2 and 
Zq is an overall constant factor which arose from the integration of the ghost fields included in the action to ensure 
that the result is gauge invariant. As we are dealing with Abelian gauge fields, that integral decouples from the others. 
From now on we omit Zq in the forthcoming calculations. C ab is the closed time path metric which reads 



1G 



C n = l 
C22 = —1 
C12 = C'21 = 

The generating function for connected n-point functions is given by 

W = -i1nZ[J,K] 

From it we define the mean values 



u SW - u sw 



and two point functions 



n (x) 4>f b (x') + G-Z (x, x') = 2C ar C bs —- 



sw 



5K l A nm (x,x'^ 



SW 



A, a (x) A vb (x') + D^ ab (x, x 1 ) = 2C ar C bs 

oKf s (x, 



x') 



& (x) A, b (x') + C W) A, b (x) + HZ (*, x') = 2C ar C bs 



SW 



(77) 
(78) 

(79) 
(80) 

(81) 
(82) 
(83) 



B. 2PI Effective Action 

The closed time path 2PI effective action is defined as 

r = W - C ab [ d^x^—g^ (x) Jl n (x) - C ab [ d^xJ—gA^ (x) J» (x) 



xK™ b (x, x') [C (x) 4>Jl (*') + G™s (x, x') 

_}_ c ar c bs J d ± x J— g J d ^ x '^gl 

xK^ b (x, X ) Afir (x) A vs (x ) + D^ivrs {x, X ) 



(84) 



where we have neglected the terms with mixed mean values, because as we shall work with a system for which the 
mean values of the fields vanish, these terms will vanish too for they couple through the mean values. 
From this action we define the sources 



ST 



<% (x) 



-C ab J—g,r b n (x) 

- l -C ar C bs J atx'^—g^' [K™ (x, x') + (x, x')] $f s (x') 



(85) 



J dV [K% (x, x') + KZ (x, x')\ A, s (x 1 ) 



ST 



SA^ a (x) 



= -C ab V~gJ^ (x) 



__ C ar C bs J d 4 x , 



(86) 



V=F [Kg (x, x 1 ) + K% (x, x')] A vs (x') 

_l_ c ar cbs J d ^^gl {Xj xl) + K ,n x , )} fa {x>) 
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sr 



5GT m (x,x>) 
ST 



c ar c bs V=gV-9'K% (x,x') 



(87) 



(88) 



Inverting these expressions and replacing them in the Effective Action and shifting the fields by 

(x) = A^r (x) + a^r (x) 



(90) 



we obtain 



r = C ap C sq / d 4 



ST 



SGpq (x,x' 



-G^" m (x,x') + 



i In / LV-Da,, exp i - iC ap / d 4 x 

JCTP J 



ST 



ST 



iC ap C sq / d 



^X / dV 



5r 



<5(/)j," (x) 
(x')^{x) + 



SA^ p (x) 
ST 



ft (x) 



(91) 



3D fii/pq (x, X ) 



(*)< (*') 



This expression can be interpreted as an implicit equation whose solution is the Effective Action. Formally this 
solution can be written as |23| 



r = S$ -ilndet [GJ" 1 ] - ilndet [G^ 1 ] - ilndet [D~l] + T 2 



(92) 



where T 2 is —ih times all the two particle irreducible vacuum graphs with lines given by Gi a b and vertices given by a 
shifted action Si n t given by 



(j) + lp 


- s 




° J dX 5K 















(93) 



C ab c a'b' / d 4 x / d i x , 



5 2 S 



i>a (x) 5<j) a > (X 1 ) 



(fib (x) (fib* (x') 



Writing 



QlbM = 



S 2 S F 



[ x )^nb i x ') 

S 2 S F 
SA^ (x) SA vh (x 1 ) ' 



(94) 
(95) 



performing the scalings 



Gib - jfG^i 

D?Ax,x ) ^^Df v {x,x') (96) 
Gtb( x , x ') ^NG lab (x,x')Sv 
GZb( x ' x ') ^Ng 2ab (x,x')S^ 
V%(x,x') -^NV%(x,x>) 

and taking into account that each loop of scalar field counts N, that each trace over scalar field indices counts N and 
that there is an overall factor of N in the action, the effective action reads 
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r = -AHndet 



(Gi)- 1 



+ -Amdct 



(G2)- 1 



+ -lndet 



+ ~Y J d 4 xy/^g J d 4 x' ^/^g'Qi ab (x,x') G lab {x,x') 
+ — / d 4 x^g / d 4 x' yf^g'Qiab {x, x 1 ) G 2a b (x, x') 



+ y J d 4 x^—g J d 4 x'^V^(x,x')Df v (x,x') 

_ N ^ C abcd J d * x ^ g ^ ix)D ab {Xi x) [Glcd ( X) x) + G2cd {Xi x)] 

+ lN ^C abc C a ' b ' c ' [ d 4 xj—g [ d 4 x'^D a ^ (x,x')g^(x)g^(x') 



(97) 



G 1W (x,x') dad a G 2cc < (x,x') 



with G ab the propagators for the scalar field, and the propagators for the electromagnetic field. 

It is easy to check that no other loops contribute to the 2PI effective action to this order in 1/N. Indeed consider a 
graph that is order N°, adding a new scalar line that maintains the 2PI condition amounts to add two vertices (and 
hence to multiply by N~ 2 ), one propagator (another A -1 ) and one momentum integration (a factor of N), which 
gives a graph that is order N~ 2 . Adding a new photon line also amounts to add two vertices, one propagator and 
one loop integration, so the new graph is again order N~ 2 . As described in Section III, G21 (x,x') is the positive 
frequency propagator and G12 (x,x') the negative frequency one. Their spatial Fourier transforms read 



G 2 i(x,x') = J 

G12 {x, x') = f d 1 /9 u* K (t) u K (t 1 ) cxp [in. (f - r')] 
J (2irr 



7~^ Uk (*) u « exp ( r ~ r ')l 
(2?r) 

d 3 K 



(98) 
(99) 



where u K (t) are positive frequency modes, that are solutions to the Klein Gordon equation for the scalar field. 

In the next subsections we shall find the evolution equations for the different two point functions involved in our 
calculations. 



C. Equation for the scalar field propagator 

The equation for the scalar field propagator is obtained by taking the functional derivative of the effective action 
with respect to the desired two point function, i.e. 



ST 



5G 



lab 



-NG-J + l -NG lab (x, x') - (x) C abcd D;l (x, x') S (x x') 



6 C acd C bc ' d 'Dli (x, x') dJ p G 2dd , (x, x') = 



2 ~ ~ » 
from where, to leading order in 1/N we have 

iG ab (x,x') = Q lab (x,x') 

Replacing the expression for Q\ ab (x,x') we obtain 

-^d^ {y/=ggT) dv + g^d^ - (to 2 + X K) ) C ab G bc (x, x') 



-J c J(x-x') 



(100) 



(101) 



(102) 



where the delta function is a coordinate function. We see that to this order in 1/N, the equations for the scalar and 
electromagnetic potential decouple and therefore there is no back-reaction of the electromagnetic field over its sources. 
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D. Equation for the electromagnetic field propagator 



It is given by 



ST 



iND^- 1 + \NVZ (x, x>) - N^- [G lab (x, x) + G 2ab (x, x)] 

2 

+iN e -C acd C bc ' d V" {x)g v0 (x') \g 1cc - (x,x')dJ G 2dd , (x,x') 



(103) 



= 



where the factor N in the term ND^ 1 comes from the scaling of the inverse of the propagator. This is already the 
equation to leading order in the 1/N expansion. Replacing T>^ (x, x'): 



Kb x ') = V^FIx) {<T (x) + (l - ±) g» a (x) g"P (x) d a d 

- -jL=d a W~gg^ a g v!i \ ApF^ - [V=g~g» a g vf3 ] A a d v A } C ab 

The equation for the propagators for the electromagnetic field then reads 

A(£ (x, x 1 ) D% (x, x 1 ) - e V (x) [G lab (x, x) + G 2ab (x, x)] D% (x, x') 

+ie 2 C acd C hc ' d ' J dx"^J-g{x")g» a (x) g" (x") 
x G lcd (x, x") dJ G 2dd , (x, x")] D% (x", x 1 ) 



S(x-x') 



(104) 



(105) 



E. Changing to conformal time 

By changing to conformal time drj — dt/a (t), using r = Hrj and performing the scalings 

G ( i\ Gj ab (x, x ) 
iab (X, X ) > 

a (t) a (t'J 

it can be checked by a straightforward calculation that the equation for the propagators read 

fi(r) 



and 



where 



and 



u * - Jp a2 W - (6X - 1) \ C ab G hc (x, x') = i5 c J (x - x 1 ) 



C ab D b J p (x,x') 



-e 2 I^ (x, x') D% (x, x>) + le ^ d C bc ' d ' J dx"YZ Ad , (x, x") D b J p (x", x') 
= iS^pSa/S (x - x') 

□x = V a0 d a dp = -d 2 a + V 2 
TZ (x, x) = rf v [G lc d (x, x) + G 2cd (x, x)} 



K V c>,dd> = vr{x)rf p {ar) [Cue {x,at>) dJ' p G 2dd , (x,x") 
To avoid over- not at ion, we have also used x as the dimensionless variables. 



(106) 



(107) 



(108) 

(109) 

(110) 
(111) 
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F. Equations for the retarded and Hadamard two point functions of the electromagnetic four potential 

In this section we shall hnd the equations for the retarded and Hadamard propagators. The latter encodes all 
the information about the state of the electromagnetic field while the former evolves the initial conditions. For the 
electromagnetic field we have 



1 



ie 2 J dx" Eff ; n (x, x" ) (x" , x') 

ie 2 J dx" E^ 12 (x, x" ) D 2 ^ (x" ,x') = iS»5 (x - x') 



(112) 



; - - ) d»d v - e 2 r^ (x, x) 



+ie 2 J dx" Eff ; n (x, x" ) D* 2 (x" , x') 
-ie 2 { rfx"E^ 12 (x, x") Dl 2 (x", x') = 



rT □* + (l - ^ - e 2 r^ (x, x) 
+ie 2 J dx" ; 21 (x, x" ) (x" , x') 
-ze 2 y dx" E^ 22 (x, x" ) C 2 i (x" , x') = 



(113) 



(114) 



^□ x + (i--)^- e 2 r^(x,x) 



+ze 2 y dx" S^ 21 (x, x" ) ^i 2 (x" , x') 

-ie 2 J dx" £££ 22 (x, x" ) L> 22 (x" , x') = - iJ^tf (x - x') 



(115) 



get 



By substracting the first two (or the last two) equation and defining £{? et = E^ 111 — S^ 212 = E 2121 — E 22 22 , we 



rf v U x + ( 1 - t ) d »d v - e 2 !^ (x,x) 



D r v e J(x,x') 



+ie 2 j dx" (x, x" ) £>"' (x" , x') = - 5£5 (x - x') 
and by adding the homogeneous equations 

D\ V1 (x, x') 

+ie 2 J dx" ££ t (x, x" ) (x" ,x') = - y J dx" (x, x" ) D% v (x" , x') 



(116) 



(117) 



wneic Zjj^ _ Zj 12 .12 ' ^21,21 — ^11.11 ' ^22.: 
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G. Fourier transformed equations 

In this subsection we find the expressions for the spatially Fourier transformed quantities that enter the Schwingcr 
Dyson equations, and for the equations themselves, in terms of the dimcnsionlcss variables defined in the body of the 
paper. The starting point is the Fourier transformation of the two point functions, i.e. 



G cd (x,x r ) = 



d 3 p 



(2tt 



,3/2 



Gcd (p, t, t') exp [-ip. (y - y')] 



where in terms of the modes the functions G c d (p, t, t') read 

G 21 (p,r,r') =/ p (r)/ p *(r') 
G 12 (p,r,r') =/ p (r')/;(r) 

For the non local kernel we have 



Xy et {X, x") 



d 3 p f d 3 q ^ 



(2^f 2 J (2tt) 3/2 ret 



S^(p,g,r,r")exp[*(p + g -).(y-r)] 



with 



(118) 



(119) 
(120) 



(121) 



£™ (p,g,T,T») - v m V ° {G\r ( P ,t,t") 8 8" G 2 n (g,r,r") - G\ 2 (p,r,r") 8 8" G 2 12 (p,r,r")} (122) 
£* et (p, g, r, r" ) = -i ( Qj - Pj ) rf°r, ij { ^ii (p, r, r" ) (g, r, r" ) - G} 2 (p, r, r" ) 5 G 2 2 (g, r, r" ) } (123) 

S;° et (p, g, r, r") = i - Pj ) r? V {Gh (p, r, r") ^oG^ (g, r, r") - G} 2 (p, r, r") d" G 2 2 (g, r, r")} (124) 



(p, g, r, r" ) = {gV + pV - pV - pV} 

x { Gh (p, r, r" ) G?i (g, r, r" ) - G} 2 (p, r, r" ) G 2 2 (p, r, r" ) } 



For the electromagnetic two point functions we write 

d 3 k 



/d s k 
- — ^D V1 (k, t, r') exp [ik. (y - y')] 



so the spatially Fourier transformed equations read 

•T - V 2 ) - (l - ±) ^ + e 2 j r^rff (p, r, r) 
-ze 2 y dr" E(£ (p, fc - p, r, r" ) (fc, r" , r') = 5^5 (r - r') 
and by adding the homogeneous equations 



D[% (k,T,r') 



D\ V1 (k,r, t') 



rT (d 2 V 2 ) " (l " J) ^ + ^ / ^fe r n (P. 
-ze 2 ! dT»XZ(p,k-p,T,r")D lul (k,r",T r ) = ^- J dr"Ef (p, fc - p, r, r") (fc, r', r") 



(125) 



(126) 



(127) 



(128) 
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VIII. APPENDIX 2: COMPUTING THE KERNELS 



A. Ultraviolet sector 



The contribution from short wavelength modes to the local term is 



d 3 p 



{2 1/3/ \f p (r)| 2 + a p p* p f 2 p (r) + &e$/t a (r)} 



(129) 



'A (2tt) 3/2 

where A is the minimum wavenumber for which the approximation of short wavelengths becomes valid. Replacing 
eqs. ([IT]) and ( |l3| ) and performing the momentum integrals we get. 

9tt 3 



If^(r,r) 



32 (2tt) 3/2 A 2 2 (2tt) 3/2 



Ci [2Ar] 



For the mixed non local: 



(p,fc-p,r,r") = 



and for the pure spatial non local 



Et'if } (P, fc - p, r, r" ) ^ 6 (r - r") 4pV 
{2|/3 p | 2 [/ 2 (r)/; 2 (r")-/; 2 (r)/ 2 (r")] 
+2a p/ 3; [/ p (r) /; (r») - f* p (r) / p (r")] f p (r) / p (r") 
2/3 P «; [f P (r) / p * (r") - £ (r) / p (r")] £ (r) £ (r")} 
Replacing the modes and Bogolyubov coefficients we are left with 

^reP (P,k-P,T,T") ~ 

9i 1 



(130) 

(131) 
(132) 



(133) 



49 (r - T")p l p l \ - — ^ sin [2p (r - r" )] + ^4 cos [2pr] - ^ cos [2 P t^ 



To evaluate the momentum integral, we must take into account that we shall be interested in the transverse component 
of the electromagnetic two point function, which is obtained by taking the curl of the corresponding equations. If we 
assume that the magnetic held propagates along the ^-direction, then the curl picks out the components x — x, x — y 
and y ~ y of the equation, and of these only the x — x and y — y are non vanishing. For either of them we get 

d 3 p 



A (2tt) 3/2 



E ret (P) (Pi k ~P, T , T " ) = 



.griaPACr-T")] 9j 

16 A 8 v n 



(134) 



3 3 

— cos f2Arl + cos f2Ar"" 

4r 4t" 



We see that for both kernels the only momentum contribution is due to the lower limit of integration and these 
terms must cancel against the ones that come from the upper limit of the integral in the intermediate momenta. 

For the pure spatial Hadamard kernel we keep only the terms quadratic in the Bogolyubov coefficients, as they will 
give the main contribution in this momentum interval, i.e. we have 



Sf (p) (p,fc-p,r,r")~4pV 
{2 |/3 p | 2 [2 \f p (r)| 2 \f p (r")| + /; 2 (r") f p (r) + f p (r") ff (r) 

+2/V* P [\f P (r)| 2 /; (r") + \f p (r")| 2 f p (t) 

+ 2p;a p [\f p (r")| 2 f P (t) + \f p (r)| 2 f p (r")] } 



(135) 



and for which the same reasoning as for the retarded kernel applies. We can conclude that the ultraviolet sector of 

the created particles spectrum makes no contribution to the noise kernels. A way to understand it is to observe that 

(2) 

the mode functions for inflation and reheating are given by the same Bessel function, Hy 2 ( z )> but with different 
arguments. This fact makes very small the amount of created particles, as can be seen from the expression for (3 p . 
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B. Infrared sector 



The local kernel for this portion of the momenta spectrum was computed in the body of the paper. In this part of 
the appendix we shall evaluate the retarded non local kernel and show that it is negligible with respect to the local 
one. The computation of the non local Hadamard kernel Ej* is straightforward. 

For the i — components of the retarded kernel we have: 



(p, k-p,T,r ,, ) = - i@ (r - t") (k l - 2p l ) 



(P* p a p - P* k _ p a k . p ) f (r) / (r) [/* (r")/ (t») - /* (r") / (r") 
+ (a* k _ p (3 k „ p - a* p (3 p ) /* (r) /* (r) [/ (r") /* (r") - / (r")/* (r") 
+ 2 (|/3 fe _ p | 2 - |/? p | 2 ) /' (r) / (r) [/ (r") /* (r") - / (r") /* (r»)] } 

while for the pure spatial components of the same kernel we get 

Ef et 0, fe - p,r,r") = 6 (r - t") (4p l p l + k l k l - 2p i k l - 2k l p l ) 

{ (|/3 P | 2 + |/3^ p | 2 ) [/ 2 (r) /* 2 (r") - /* 2 (r) / 2 (r")] 

+ (ajk-p^^ + ap^) [/ (r) /* (r") - /* (r) / (r")] / (r) / (r") 
+ (/3 fe _ p a^_ p + [/ (r) /* (r") - /* (r) / (r")] /* (r) /* (r") 



(136) 



(137) 



Recalling again that we shall take the curl of the equations and considering for example, that the magnetic field 
propagates along the z direction, we need only the — x and — y components of the mixed kernel and the x — x, 
x — y and y — y of the pure spatial one. But except for the x — x and y — y of the pure spatial noise kernel, the rest 
of the components, by virtue of the momentum independence of the mode function (|18|) , are odd functions of p and 
hence their momentum integral vanishes. We then find the desired decoupling of the equation for the transverse part 
of the pure spatial propagator. Of course the equations for longitudinal part of the propagators are still coupled and 
carry the information about charge conservation. 



Performing the momentum integral of the component E,. ei 

r A _ d 3 P 

Jt 



x(N) 



(2tt 



,3/2 



for example, we get 

£^ P) (p,fc-p J T J T*)~Tre(T- 



8 A 2 

3 A 



(138) 



2 3/7-4p2 



3 + 2 7 
27 



z 



H 



3/2 7 ( Z ) ^3/2 7 



5/7 
z 

H?L (z" 



3/2 7 



'#§27 K 



H. 



3/2 7 ( Z ) ^3/2 7 ( Z ") 



H. 



3/27 ( Z ")^3/27 ( Z ) 



(2) 
3/2 7 



(1) 



r(i) 



3/2 7 



3/2 7 



3/2 7 



Comparing the prefactor in this equation with the corresponding one in eq. ( p3| ) we see that by virtue of the 
logarithmic dependence in the latter, the non local kernel turns out to be several orders of magnitude smaller than 
the local one. 
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